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Abstract guidelines of the Investigation, and some of Its 


Effects of three different models for the 
treatment of subsonic boundary conditions, applied 
to the problem of flow in a channel with a bump, 
ore discussed. A preliminary discussion of the nu- 
merical treatment of the corners is presented. 

1. Introduction 

Experimentalists have learnt to appreciate 
the importance of a proper installation, a careful 
calibration of instruments and an analysis of en- 
vironmental interferences. „ Data obtained by the 
Wright Brothers on their historical wind tunnel en- 
abled them to evaluate the possibility of flying; 
by current standards they would probably not even 
fit in what engineers mention as "ballpark" range. 
The designer's problem is efficiency; errors larger 
than IS in basic information should not be tolerat- 
ed. 

Surprisingly, a large amount of numerical 
work is far from reaching such standards. More 
surprisingly, very little seems to be done to im- 
prove the situation, and this is so much more dis- 
turbing since some claims are occasionally made 
that nume-ical analysis should replace experimei.ts. 

The present paper is a modest attempt to 
show a possible approach to understand reliability 
of numerical analysis. It is shaped as a series of 
numerical experiments. Empiricism, however, is not 
suggested. Empiricism is manifested by juggling of 
arbitrary coefficients, mesh refinements and, mani- 
pulation of arbitrary additional terms. The work 
described in this paper, instead, is inspired by 
the idea that a numerical procedure describes, more 
or less accurately, a physical model and that the 
understanding of such a model will lead us to Judge 
whether or not our calculation makes the physical 
sense which it should. To go back to our former 
comparison with experiments, a typical wind tunnel 
correction inspired by the same criterion to which 
the present study is inspired is the wall correc- 
tion for viscous and transonic effects. 

In preparing the present paper, a very 
large amount of cases were computed, expressing 
different lines of tnought and, for each one of 
them, tests of different parameters. A detailed 
discussion transcends the limits of a presentation. 
Therefore, I will limit myself to showing the 


highlights, without attempting to be exhaustive and 
even to draw conclusions which, as it will appear 
from the context, could, at this stage, still be 
hasty and Inappropriate. 

<?. A channel flow 

On September 18-19. 1979, a Workshop was 
held in Stockholm, the object of which was the com- 
parison of results obtained by using different nu- 
merical methods on two assigned problems, the 
second of which was formulated as follows. 

"Internal two-dimensional flow through a 
parallel channel having a 4.2S thick circular arc 
•bump' on the lower wall. The ratio of static 
downstream pressure to total upstream pressure is 
0.62351 2 (corresponding to M=0.85 in iser.tropic 
flow), and the distance between the walls is 2.073 
times the chord length of the bump." 

Obviously, the emphasis of the assignment 
was on steady solutions and transonic flow with an 
imbedded shock. The latter requirement adds a 
number of complications to the problem of a subson- 
ic, steady, isentropic flow in a channel. The as- 
signed data were so close to producing a choked 
flow that some of the methods generated a choked 
flow (all potential fully conservative methods) and 
others did not (all potential nonconservative 
methods and Euler solvers, and Hafez's artificial 
compressibility method). Scattering of results and 
conflicts between conclusions are not new in our 
short history of numerical analysis. As I recall, 
the first numeric?! contest was inspired by Morton 
Cooper in 1965, a calculation of blunt body shock 
layers for an ellipsoid of revolution with a 2: 1 
axis ratio, at a free stream Mach number of 3 [1]. 
Techniques ranging from truncated series expansions 
to integral relations to inverse methods offered a 
variety of results. Comparing them with what has 
now been accepted as standard, that is, a second- 
order finite difference calculation with bow-shock 
fitting, we see that methods focussed on the stag- 
nation line gave good results near the stagnation 
line and poor results away from it, whereas methods 
focussed on the sonic line had the opposite 
behavior [2). The contest clea* ly showed a need 
for a different numerical approach, more general 
and powerful. 

The object of a Workshop is, indeed, to 


promote healthy competition and unrestrained dc- 
ba< •csj not to solve problems or to emit verdicts; a 
Workshop can be considered successful If It In- 
spires new, and deeper, work. In studying the 
results of the Stockholm workshop, I decided to 
take a closer look at the channel flow, at least 
Unitedly to a certain brand of Euler solvers. 


where values at three adjacent points, from the 
boundary In, are denoted by 1,2 and 3, sequential- 
ly. 

J3. Computational grid 


Obvious questions to be answered were: 

1) Is a steady state reached? 

2) Do results depend on the type and size 
of the computational mesh? 

3) How do different treatments of the left 

and right •' boundary of the computed region affect 
the results? . 

4 ) Can any detail be provided of the flow 
near the leading and trailing edges of the bump? 

Kote that I abstain^ from mentioning in- 
tegration schemes. Relative virtues and shortcom- 
ings of such schemes, including their ability to 
capture shocks, their numerical diffusion and 
dispersion, etc. are out of context. They cannot 
be tested In the channel flow problem unless the 
questions above have been exhaustively answered. 
On the other hand, there are general features of 
the flow which should be revealed and which should 
provide clues to the questions, regardless of the 
Integration scheme having been used, at least so 
long as the flow 1* far from transonic. 


The grid suggested for the Stockholm 
workshop was a Cartesian grid, normalized between 
upper and lower wall of the channel, and stretched 
from — to with an strong accumulation of grid 
lines over the bump. I used this type of grid, 
forcing two grid lines to originate exactly at the 
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leading and trailing edges (Fig. 1). 

I also adopted a different grid (shown In 
Fig. 2) which Is obtained using a conformal mapping 
of the Karman-Trcfftz type: 


i I decided, thus, to limit a preliminary in- 
vestigation to subsonic, isentroplc flows, ar.d I 
adopted the MacCormack, predictor-corrector scheme 
to the equations of motion In the form: 
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where z=x+ly Is the complex coordinate In the phy- 
sical plane, c = (*in is the complex coordinate In 
the mapped plane, and 4 is related to the 
thlcl:ness-to-chord ratio of the bump, i, by: 


V *(V.V)V ♦ T VP = 0 
t 

* 

where P Is the logarithm of the pressure, V the 
velocity, T the temperature and y the ratio of 
specific heats. Pressure and temperature are re- 
lated by 


P S - 1 - In T (2) 
Y-1 

For subsonic flows, the MacCormack scheme Is safely 
applicable to (1) and has the advantage of great 
simplicity. To maintain second-order accuracy at 
the boundaries, where the MacCormack 3chemo can be 
applied only at the predictor (or corrector) level, 
for want of external data in the other level, any 
derivative at such a level is discretized by 
differences of the type: 
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The n-coordlnates arc normalized between the upper 
and lower walls, and the (-coordinates are 
stretched in the horizontal direction a3 In 'he 
Cartesian grid. Calling u and v the velocity com- 
ponents in the direction of the Cartesian axes for 





the first grid, and In the direction of the 
t = constant, n » constant axes for the second 
grid. It turns out that v vanishes exactly along 
the upper wall when the first grid Is used, and 
along the lower wall, when the second grid is used. 
Consequently, one may expect the calculation to be 
easier and, perhaps, the results to be better along 
the wall where v vanishes. In any event, the boun- 
dary conditions at rigid walls arc enforced by 
first Integrating the Euler equations as at Interi- 
or point- , .nd then by correcting the pressure to 
satisfy r ie vanishing of the velocity component 
normal to the wall [ 3 , ** 3 . Along walls wnere v Is 
not Identically zero, the u-momentum equation Is 
replaced by an equation along the tangent to the 
wall. 


Internal flows, several simple physical models can 
be adopted, of which here Is a sample: 

1) The region of Interest Is a channel of 
a finite length, connecting two Infinite cavities; 
the gas Is at rest everywhere; at t=0, the stagna- 
tion pressure is increased In the cavity at left, 
until a given value Is reached, and then kept con- 
stant, 

2) The same setting Is used, but at t=0 
the pressure In the c ivlty at right Is decreased 
until a given value Is reached, and then kept con- 
stant, 

3) The channel Is Infinitely long and It 
contains- a gas at rest; at tsO, the channel Is ac- 
celerated towards the left, until a cruising speed 
Is reached. 



Calculations are actually performed In a 
computational plane, (X, Y), where the grid Is 
evenly spaced In both directions. If the first 
grid Is used, additional coefficients appear, con- 
taining dX/dx, 8 Y/3x and 3Y/3y. If tnc second grid 
Is used, two Independent sets of ••dditional coeffi- 
cients and terms appear, the first due to normali- 
zation and stretching, and containing 
dX/dc, 3Y/3£, 3Y/3n, the second due to the mapping, 
and containing 
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The major physical difficulty Is offered by 
the leading and trailing edge corners, where the 
flow stagnates. In using the first grid, the dif- 
ficulty Is reflected In the discontinuity In the 
slope of the lower wall (which affects all points 
on the vertical grid lines Issuing from the 
corners). If the second grid is used, a mapping 
singularity appears at the corners and the equa- 
tions of motion, expressed in terms of C and n, be- 
come indeterminate. In both cases, thus, some spe- 
cial treatment must be given to the corner points 
and their Immediate neighbors. If the grids are 
laid to avoid passing through the corners, the ef- 
fect of neglecting them has to be evaluated. 


5.* Inlet and outlet boundary conditions 

Another critical Issue regards the treat- 
ment of the arbitrary computational boundaries 
which delimit an Inlet and an outlet to the region 
of interest. Such boundaries cross regions of sub- 
sonic flow and some physical model is required to 
supply the Information from outside which Is neces- 
sary. New Interest seems to have arisen on this 
problem In recent times, but the physical Implica- 
tions of model lrg a subsonic boundary seem not to 
have been grasped firmly yet. The problem of sub- 
sonic boundaries cannot be disassociated from the 
problem of choosing Initial conditions [5). In 


In the first two cases, two models of tran- 
sitions from the interior to Infinite cavities are 
adopted at each boundary point on the left and on 
the right. As explained In [5], one can stipulate 
that the fictitious flows in the transitions are 
quasi-steady (the length of the transition being 
assumed as vanishingly small), so that, for the 
purpose of closing the boundary data sets, steady 
equations of motion can be differentiated in time. 
G.i the left, the total pressure and the slope of 
the velocity vector at each entry point are as- 
signed. The latter condition brings In the largest 
arbitrariness In the model. Physically, one can 
always Justify a choice of slopes by assuming that 
the Inlet Is equipped with a series of guiding 
vanes. In the present case, for example, one can 
assume that all velocity vectors are parallel to 
the rigid walls; this is obviously not the case for 
an infinitely long channel, and the effect of such 
a restrictive assumption on the rest of the flow 
has to be evaluated. 

* 

The equations used at the Inlet are: 

1) the definition of total pressure, dif- 
ferentiated In time under the assumption that the 
total pressure itself may be a function of time: 

2 

TP ♦ u (1 ♦ o ) u » T P (6) 

t t o ot 

Where the Index o denotes stagnation values In the 
Infinite capacity, and o - v/u is a prescribed 
value , 

2) a left-running characteristic equation: 

a P - t u = It (7) 

V v 

where n Is the left-hand side of (7) as computed by 
the standard integration routine. 

The outlet model Is simpler, since the v- 
component of the velocity Is determined on the 
basis of Internal Information only [OJ; In tho 
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present case, It Is sufficient to prescribe the 
•lit pressure as the pressure In Infinite capacity, 
and compute the u-compor.cnt of the velocity accord- 
ingly. The equations are, thus: 

1) the contlnutlty equation: 


T 


PUP 4 TP U spuP ♦ TO U 

t t • - -t « »t 

2) the definition of total pressure: 
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3) a right-running characteristic equa- 
tion: 


a P ♦ y u = R (10) 
t t> 

where R has the same meaning as In (7). In (8), 
(9) and (10), v Is computed by the standard 
routine: u is unknown but It can be eliminated 
easily. Naturally, here too there is an element of 
arbitrariness, whose effects have to be checked. 
For example, in an infinitely long channel, the 
pressure across the channel Is not exactly constant 
at a finite distance from the bump. 

The Inlet and outlet boundaries Just 
described allow perturbations proceeding from the 
Interior to interact with the conditions In the in- 
finite cavities. For each perturbation reaching 
the boundary, a new perturbation Is generated and 
transmitted In the opposite direction. The process 
will eventually reach an asymptotic steady state, 
but the number of waves of a sizeable amplitude 
moving back and forti: can be very large. 


The third model relies on a sluple Idea: if 
the motion were one-dimenslor.al , all perturbations 
would travel outwards as simple waves, at the end 
of the acceleration phase. A simple wave Is easy 
to describe using information from the interior and 
the constancy of one Riemann invariant from the ex- 
terior. In a two-dimensional problem of internal 
flow, the wave3 cannot be exactly simple waves, but 
no major errors are expected if the velocity vector 
is forced to be parallel to the rigid walls at the 
inlet. The simple wave equations at the inlet are 
modified as follows: 


a P - t (1 «■ o^J^^u = R 
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Similarly, at the outlet: 

aP ♦ ru u /q = R 
t t 

( 12 ) 


■^t “ YU u fc /q * tv v fc /q 

where v^ Is computed by standard routines and q Is 
the modulus of the velocity. 

5 . Two-dimensional calculations with models 1_ and 2 

We describe now the general features of 
calculations made using the first two models men- 
tioned in Section l '. One of the problems presented 
by a study like the present one is the large amount 
of data produced by a single run and the necessity 
for organizing them in a series of simple plots, 
easy to Interprc" . 1 decided to store the follow- 
ing Information: 

1) at every step, P and u on the lower 
wall, at the Inlet, at two selected points, and at 
the highest point on the bump, 

2) at every step, location of selected 
isobars on the lower wall (to build an Isobar pat- 
tern on an (x.t)-piane) , and > 

3) at selected steps, P, u and v at all 
the grid points; this Information can be easily 
processed to provide Mach numbers and total pres- 
sures. 

The basic geometry has been defined as a 
channel with a width equal to 2, containing a tump 
which extends from x=-1 to x=1 and which has a max- 
imum thickness of 0.2. This defines a corner angle 
of 157 . 38 °, and 4 = 0 . 87 ^ 3 . To avoid initial compli- 
cations at the corners, so that our attention can 
be focussed on the wave propagation and the effects 
of boundaries, we use a smooth lower wall which can 
be easily obtained from the mapping function by de- 
fining the wall as the image, in the z-planc, of a 
line n * b, where b is a constant greater than 0. 
The same definition can be transferred to the code 
which uses the Cartesian grid. In Fig. 3 there are 
some shapes of the lower wall for different values 
of b; one can observe that, for b less than 0.01, 
there is no practical difference between the wall 



Fig. 3 

so defined and the wall corresponding to b=0. 


Vc begin with a very smooth wall, defined 
by b*0.1. The computational mesh has 7 Intervals 
between the rigid walls and 30 Intervals In the x- 
• direction, stretched between x*-2.3A5 and x=2.3A5; 
16 Intervals cover the bump region. The stagnation 
pressure Is raised (In the first method) or the 
exit pressure Is lowered (In the second method) to 
produce final values of the Mach number "at Infini- 
ty" of the order of 0.1. A plot of P vs. time at 
the 6th node on the lower wall Is shown In Fig. A, 
for the ease where the first mchod Is applied. The 
oscillations arc obviously produced by waves going 
back and forth along the channel (at such low Mach 
numbers, the spied of propagation Is practically 
the same In both directions, and the phenomenon 
. shows a well defined frequency); one disturbing 
feature of the model Is the smallness of the damp- 
ing factor. 



Fig. 6 

A similar pattern appears (Fig. 5) using of level lines at a given step. For example, lso- 

the second model; the oscillations are smaller In bars and lsomachs at step 1000 (Figs. 7 and 6) look 

amplitude but still damped very slowly; Fig. 6 very reasonable, although the Isobars would r.ot 

shows u(t) at the same node ano from It we see that pass a closer scrutiny, due to a clear lack of sym- 

thc velocity presents smaller oscillations but that metry. As a matter of fact, If we plot T at the 

o' steady state Is far from having been reached 6th node vs. P at the ?MLh node as they evolve In 

after 2000 computational steps (t=1^6). time, we sec that, after 2000 steps, the plotting 

line still oscillates between -0.00283 and 
-0.00353. whereas at both points P should be about 
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Such details ere hard to detect from plots 
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Fig. 8 


-0.0032. 

An analysis of these oscillations should 
take at least two elements Into consideration: the 
first Is the Mach number effect, and the second is 
the influence of geose-ry. To nave an idea of the 
Mach number effect, let us rerun the above cases 
•for a Mach number of 0.5. Plots of P and u vs. 
time arc shown in Figs. 9 and 10 for the first 
model; Fig. 9 should be compared with Fig. 4. Os- 
cillations still appear but they seem to be damped 
much more quickly. A similar behavior is seen in 
Figs. 11 and 12, which refer to the second model. 
One should, however, take care not to draw hasty 
conclusions from Figs. 4 and 9. or Figs. 5 and 11. 
The scale of P in Fig. 9 is 20 tines smaller than 
in Fig. 4, and in Fig. 11 is 40 times smaller than 
in Fig. 5; the correct conclusion is that very 
small pressure waves take a long time to be elira- 
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inated. The patterns of level lines (Isobars in 
Fig. 13, isomachs in Fig. 14) are much better than 
their counterparts for M=0. 1 (Figs. 7 and 8). Even 
the v=constant lines, which are very critical, look 
good (Fig. 15). At this stage, it pays to take a 
look at lines of constant stagnation pressure (Fig. 
10); here a new element appears. In fact, the 
Stagnation prt.jute is practically constant every- 
where, but it drifts away in the vicinity of the 
"corners" (or whatever remains of them in the 
smoothed wall). The stagnation pressure Is a very 


lens'*' parameter, Indeed, and It Is the proper 
Indicator of local inaccuracies, when a steady 
state Is apparently reached numerically. In this 
ease, It Is obvious that Inaccuracies should be at- 
tributed to the vicinity of a singular point of the 
napping and to the consequent worsening of the 
metric. 
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Before going Into the difficult problem of 
the corner singularity. It Is proper to develop 
more familiarity with the wave propagation pattern 
for models 1 and 2, and their possible relationship 
with the existence of a bump. We have seen, so 
far, that waves tend tc continue swaying back and 
forth, with very little damping, at low Mach 
numbers. To Judge whether the geometry, and par- 
ticularly the presence of a bump, has anything to 
do with the wave behavior, we examine two cases. 
One Is an obvious choice, a straight channel with 
no bumps. It can be easily obtained from either 
the code using the Cartesian grid or the code using 
the mapping by setting the thickness of the bump 
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Fig. 14 
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Fig. 18 
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equal to zero. 

Here are some results, for comparison with 
the previous cases. The first (Figs. 17 and 18) 
uses model 2, with a Mach number of 0. 1. Compare 
Fig. 17 with Fig. 6 and Fig. 18 with Fig. 6. The 
second uses model 2 again, but with a Mach number 
equal to 0.5 (Figs. 19 and 20). Compare Fig. 19 
with Fig. 11 and Fig. 20 with Fig. 12. Note that 
in this ease the steady flow In the channel Is uni- 
form, with n pressure equal to the exit pressure; 
the transition from stagnation pressure (PrO) to 
the channel pressure takes place In the fictitious 
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transition which has been modeled nt the Inlet. In 
Figs. 21 (to be compared with Fig. 9) and 22 (to be 
compared with Fig. 10) the first model Is used. 
Here the stagnation P Is raised to a positive value 
•nd the exit P remains equal to zero; the latter 
Is, thus, the asymptotic value of P in the whole 
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and that they have a physical Interpretation of 
their own, unrelated to the geometrical complica- 
tions of the channel. 

6, Two - dimensional calcul atlons with model 3. 

We expect calculations made using the third 
model of Section ** to converge to a steady state 
much faster than the previous ones, since the ini- 
tial perturbation affects the entire flow field and 
whatever is not pertinent to the final state is 
promptly eliminated through the boundaries which, 
in this aodcl, arc not reflective. The expecta- 
tions arc confirmed by Figs. 23 and 2'M (which 
should be compared with Figs. *1,5 and 6, respec- 
tively. Note also. In Figs. 25 and 26 (isobars and 
isomachs, respectively) how close the pattern is to 
the symmetric pattern of a steady state; compare 
these figures with Figs. 7 and 8. We omit present- 
ing results for M=0.5; they are equally good and 
not dissimilar from the ones obtained using the 
second model, although a close inspection may re- 
veal some advantage in using the third model (for 
example, the v=constant lines appear more symmetri- 
cal than in Fig. 15). 

At this stage, we conclude that: 

1) any one of the three models is accept- 
able as the description of a physical evolution, 

2) any one of the three models is accept- 
able for the evaluation of a subsonic, steady 
state, 

3) In any of the three models the computa- 
tional region can be limited to a small portion of 
the channel, bracketing the bump, 

*1) the third model, however, provides fas- 
ter convergence to a steady state, particularly for 
low Mach numbers. 

Therefore, all further investigations wll ) 


All these patterns are similar to the ones 
with the bump; we can conclude that the oscilla- 
tions arc produced by the models of the boundaries 


7.. The corner probl em 
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We may turn now to the analysis of the dif- 
ficulties arising when the lower wall has an actual 
corner. Probably, the first Idea which may come tof 
one's mind consists of performing successive calcu- 
lations, with all parameters unchanged except b. 
Let us recall that all the preceding exercises have 
been performed with b=0.1, which corresponds to a 
rather smooth lower wall. What happens If b Is re- 
duced to 0.01, 0.001, and so on? 

For b=0.01, and M=0. 1, the P(t) and u(t) 
patterns do not show sizeable difference from the 
ones of Figs. 23 and 2U (b=0.1); the lsooars and 
lsomachs plots, however, start showing signs of de- 
generation near the corners, particularly the one 
on the left (Figs. 27 and 28). The region of high 
pressure tends to spread to the left, ar.d the 
minimum Macn number Is definitively misplaced. 
Minor changes In the rest of the plots can be as- 
cribed to consequences of the poor accuracy near 
the corners. Although the calculation can be pci — 
formed for much smaller values of b (even for 
b=0.0005) without catastrophles, such results can- 
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Fig. 26 

conducted using the third model. 


corner Is singular, so long os the corner angle Is 
larger than 90 . Since the velocity vanishes at 
the corner, It will behave, In Its Immediate vicin- 
ity, in the velocity of an Incompressible flow. 
How, for the present case, the velocity of an In- 


compressible flow nlong the lower wall Is propor- 
i .t tonal to the G of (5). To learn more about u and 

Its dependence on b, let us plot. In Fig. 29, o 
i , aeries of curves of G(x) for different values of b. 

In the sam« figure, arrows Indicate the approximate 
, . positions of the nodal points for the computations 
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shown In the preceding Sections. 

The qualitative behavior of u Is expected 
to follow the trends of Fig. 29, even in a compres- 
sible, unsteady flow. Similar considerations can 
be developed for P. At the corner Itself, the only 
quantity' to be computed Is P, since u -jr.d v vanish 
Identically. For all practical purposes, a reason- 
able estimate of P is obtained from the definition 
of tne stagnation pressure, applied to values aver- 
aged on the 3 points surrounding the corner. Such 
an estimate tends to become more and more accurate 
as the flow approaches a steady state. The diffi- 

* 

culty, thus, does not lie at the corner but at the 
nclghborirg points. From Fig. ?9 it is clear why a 
derivative such as u^ may be approximated by a 
2-point difference if b=0.1 but the approximation 
is very poor if b=0.01 and it becomes disastrous if 
b=0. On the other hand, the values at the corner 
point cannot be Just skipped because, if the 
derivatives at the neighboring points are approxi- 
mated by onc-sldcd differences only, all connection 
between the two sides of the corner is lost and, if 
the derivatives are approximated by differences 
between the two neighboring points, again the ap- 



Fig. 30 

"proxlmation is poor, as we can see from Fig. 29. 


values of the parameter b. Again, for b=0 or close 
to it, the values at the corner can be determined 
without difficulty but the derivatives at the 
neighboring points are hard to approximate. Cot se- 
quences similar to the ones observed In preceding 
Sections arc visible in Figs. 31 and 32, where 
lsomach lines arc shown for the cases, b = 0. 1 and 
bsO.Ol, respectively. 

Fortunately, the Influence of the corner 
singularity is limited to a very small region sur- 
rounding the corner; in the grid used in our calcu- 
lations, the departure of G from the very smooth 
curve, relative to b=0.1, occurs only inside one 
cell and differences between the curves relative to 
b:0.01 and b=0 appear only in a negligible portion 
of the cell. The motion In the neighborhood of the 
corner must be considered as C.ie same which would 
take place in the presence of a fairing, such as 
the one defined by bsO.I, ,1.us a sort of triangular 
region, comprising the corner, where the flow stag- 
nates in a succession of quasi-steady states, even 
if the general motion is unsteady. Such a region 
should be considered only to evaluate the pressure 
at the corner, which should not be used, however, 
to compute derivatives at the neighboring points. 
It seems that analytic expressions for P and u 
within the corner cell, based on "reasonable" solu- 
tions for incompressible flow fall to provide the 
proper approximations to the derivatives, oerhaps 
because their domain of validity is too small as 
compared with the size of the cell. A safer ap- 
proach could consist of replacing ti * values of g 
and $ at the corners with their values at the 
corresponding point defined by b=0.1, say, in this 
way providing a localized fairing which should re- 
flect the physical behavior, maintaining a geometry 
which would not conflict with continuity require- 

C» 

ments and producing values for pressure <*nd veloci- 
ty components which could be used directly for the 
approximation of derivatives. More on this sub- 
ject, including a study of different corner angles, 
down to the limiting (and very special) case of a 
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Fig. 31 
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yu corner, will appear in a forthcoming paper. 


If the Cartesian grid is used, the dlltl- 
culty obviously persists. Plots of the slope of 
the lower wall arc shown In Fig. 30 for various 




Fig. 32 

• 

Acknowled g mcnt . 

Part of this work was sponsored by the NASA Langley 
Research Center under Grant No. NSG 1298 and per- 
formed during the Author's visit at ICASE In August 
1980. 


References . 

• 

1. Perry, J., "A comparison of solutions to some 
blunt body problems," ATAA J. JJ, 1925-1926 (1966) 

2. Morcttl, G. and Abbett, M., *-A time-dependent 
computational method for blunt body flows," AIAA J, 
A. 2136-2191 (1966) 

3. Abbett, H.J., "Doundary condition calculation 
procedures for lnvlscld supersonic flow fields," 
Proc. AIAA Comp. FI. Dyn. Conf. , Palm Springs, 
153, 1973 

A. deN’eef, T. , "Treatment of boundaries In usteady 
lnvlscld flow computations," Delft Unlv. of Tcchn., 
Dept, of Aerosp. Eng. Rep. LR-262 (1978) 

5. Korettl, G. and Pandolfl, M., "On the calcula- 
tion of subsonic flows In ducts,", POLY H/AE Rep. 
No. 80-18 (1980) 

6. Pandolfl, M. and Zannettl, L., "A physical ap- 
proach to solve numerically complicated hyperbolic 
flow problems," VII Int. Conf. on Num. Mcth. In 
Fluid Dyn. (I960) 

1 . Van Dyke, M.D. , Perturba i on Me thods In Fluid 
Meehan 1 c? , Parabolic Press (1976) 


■ i'-j 


